On the existence of internal modes of sine-Gordon kinks 
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We study whether or not sine-Gordon kinks exhibit internal modes or "quasimodes." By con- 
sidering the response of the kinks to ac forces and initial distortions, we show that neither intrinsic 
internal modes nor "quasimodes" exist in contrast to previous reports. However, we do identify a dif- 
ferent kind of internal mode bifurcating from the bottom edge of the phonon band which arises from 
the discretization of the system in the numerical simulations, thus confirming recent predictions. 
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I. INTRODUCTION 



Solitons, which were originally a concept arising in the 
study of integrable systems , have rapidly become a key 
paradigm in perturbed, non-integrable systems describ- 
ing very many physical applications ^| . A very succesful 
picture of soliton dynamics under different perturbations 
has resulted from the use of collective coordinate (CC) 
techniques In many instances, solitons behave ba- 
sically like point-like particles, and therefore their time 
evolution can be predicted by equations of motion which 
are ordinary differential equations. Nevertheless, along 
the last decade it has been progressively realized that in- 
ternal degrees of freedom of the solitons play a crucial 
role in a number of problems: They govern resonant soli- 
ton (or solitary wave) collisions Q] , give rise to nontrivial 
soliton- impurity interactions , and can be excited both 
by ac forces || and by thermal noise jjj. Therefore, it is 
very important to assess whether or not physically rele- 
vant solitons possess internal modes. 

In this Rapid Communication we study solitons of 
the sine-Gordon (sG) equation, which arises in a diverse 
range of areas of physics, covering from crystal dislo- 
cation theory to Josephson junctions flj-[|. In spite of 
the fact that linear stability analysis leads to only one 
eigenvalue in the discrete spectrum, corresponding to a 
zero frequency Goldstone mode Q , Boesch and Willis ||] 
claimed that they found an internal "quasimode" above 
but close to the lower phonon band edge. They described 
this "quasimode" as a long-lived oscillation of the width 
of the sG kink. The possibility of such a mode had been 
suggested earlier by Rice by means of a variational 
approach that reproduced with great accuracy the inter- 
nal mode frequency of the 4 model [jlj and predicted the 
existence of a similar mode for the sG problem. To our 
knowledge, Ref. Q is the only paper reporting the obser- 
vation of this mode, although hints of its existence are 
scattered all over the soliton literature . A few years 
later, it was theoretically found [pd], within the Rice ap- 



proach, that the internal mode of the sG and the (j> 4 
kinks could be excited by a constant external force or by 
an initial distortion; however, no numerical evidence was 
presented to confirm these results. Therefore, we decided 
to revisit the issue of the existence of the internal mode of 
the sG equation in order to give a definitive answer to it. 
Our approach, based in part on our earlier findings |p|, 
overcomes the difficulty that the internal "quasimode" 
would lie within the phonon band by allowing us to look 
at driving frequencies well in the gap. As we will show 
below, this permits us to establish clearly that this pu- 
tative internal mode does not exist. Furthermore, we do 
observe the excitation of a discreteness induced internal 
mode ]l^ , p^[ , whose frequency is in excellent agreement 
with recent predictions by Kevrekidis and Jones 14 1. We 
substantiate these claims by briefly recalling our main an- 
alytical results in Sec. [n], collecting our numerical results 



in Sec. Ill, and discussing our conclusions in Sec. IV 



II. RESONANCES DUE TO AC FORCES 

We begin by considering the perturbed sG equation 

<t>tt - 4>xx = - sin</> + esm(St + S ). (1) 

where e, 5, and i5o are the amplitude, frequency, and 
phase of the external periodic force. In || we ana- 
lyzed analytically and numerically the above problem in 
the more general context of nonlinear (possibly damped) 
Klein-Gordon systems, which includes Eq. (Q) as well 
as the (ft 4 model. Particularizing our results, obtained 
by means of the generalized travelling wave Ansatz fig] 
combined with the Rice Ansatz we find that the mo- 
mentum P(t) and the width of the kink l(t) (of a general 
nonlinear Klein-Gordon equation; specific constants for 
the sG case will be given below) obey the equations 



dP 
~dt 



= -qesin(St + S ), P(t) 



MqIqX 
l(t) ■ 



(2) 
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-1, 
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where X(t) is the kink position. For sG, q = 2-it, M$ = 8, 
lo = 1 and a = 7r 2 /12. The equation for the momentum 
is linear and can be solved exactly, yielding 



P(t) 



P(0) 



g e cos(i5o) 



+ V cos(<5t + (J ) (4) 



We thus see that the term P{t) 2 in Eq. (||) involves two 
frequencies, S and 2(5. As shown in for the 4> 4 kink 
this leads to dramatic resonance effects when the ac driv- 
ing frequency 5 is half the frequency £li of the internal 
mode. Without resonant excitation of kink modes, the 
kink oscillates around its initial position as a pure point- 
like particle, as shown analytically and numerically in 
p6| . However, as shown in 0] for the above parame- 
ters, in case an internal mode is excited as in the <f> 4 
system, the system energy grows and eventually leads to 
chaotic motion of the kink. In addition, the resonance 
at Cli/2 turns out to be much stronger than that at fli, 
which makes it a convenient way to probe those modes: 
In case the sG kink would possess such a degree of free- 
dom, we should observe similar phenomena. In ||, the 
frequency of the numerically observed "quasimode" was 
uj s = 1.004 ± 0.001 in our units, whereas Rice's theoret- 
ical analysis leads to Q 0^ = 1.103. In both cases, the 
corresponding half frequencies lie well within the gap; we 
can thus force the system at the predicted resonances 
without exciting much radiation which could mask any 
internal mode effects. 



III. NUMERICAL RESULTS 

In this Section, we present the results of our numer- 
ical search for internal modes or "quasimodes" of the 
sG equation. We have computed the numerical solution 
of the perturbed sG equation (|l|) by using the Strauss- 
Vazquez scheme |rit| , the total length of the system be- 
ing L = 100, the steps being Ax = 0.05, At = 0.005 (or 
other values when indicated). Free boundary conditions 
were also used, and the final time in the simulations was 
25 000, enough to detect any possible resonance M. 



A. ac force 

First, we have perturbed our system with an ac force 
of amplitude e = 0.01 and phase So = ir/2. We have 
used a sG kink at rest as initial condition. We have 
computed the energy (as defined in jlTj ) during the evo- 
lution, and then the mean value of the energy in the time 
interval 10 000 < t < 25 000, studying a range of values 
of S around half of the possible internal mode frequen- 
cies. Figure |l| collects our results for the mean energy 
as a function of the frequency. We clearly see that the 



mean energy rapidly increases for some values of 6 which 
do not coincide with any of the proposed internal modes 
or "quasimodes." The first (leftmost) maximum lies be- 
low half the frequency of the lower phonon band edge; 
therefore, this one can not be associated with a phonon 
mode and we will discuss its origin in detail below. As 
for the other peaks, Table || shows that we have been 
able to identify them with half of the values of the first 
frequencies of the phonon band, w„, which are given by 



1 + [(2tto)/L] 2 , n = 1, 2, 3, . . . , N = L/Ax. (5) 



A comparison of the third and fourth columns of Table 
H where we list the half values of the frequencies cor- 
responding to the first 8 radiational modes and the nu- 
merical peak frequencies uj n in Fig. |l|, counting from the 
second maximum from the left, allows to verify that these 
values are nearly identical in all cases (see also the rela- 
tive differences |w n /2 — w n \ in Table |). We stress that 
identification of radiation modes as the origin of these 
resonances makes perfect physical sense: Indeed, an ex- 
tended CC approach including radiation modes leads to 
the conclusion that the radiation modes should also be 
excited parametrically by ac drivings. In view of all this, 
the main conclusion we can draw from Fig. [l] is the ab- 
sence of any resonance due to the internal quasi-modc 
in the phonon band reported by Boesch and Willis [||, 
which, in view of the frequency they reported, would lead 
to an extra peak located between the two rightmost ones 
in Fig. [j]. In this regard, it is important to note that the 
accuracy of our simulations in detecting every radiation 
mode rules out the possibility that we have missed this 
quasi- mode (at either uj s or CIr). 

Interestingly, we do observe an internal mode: Let 
us recall that the first maximum of the energy occurs 
at S = uj d = 0.4999, i.e., 2u d = 0.9998 is below the 
lower phonon frequency uj p h — 1 as we mentioned above. 
The reason for this peak is a novel internal mode of the 
sG kink, that appears as a bifurcation from the lowest 
phonon of the continuum equation, of frequency ui p h = 1 , 
due to the effect of discreteness ]l2|~|T^] . As has been 
recently shown jL4j, the frequency of this new internal 
mode is uj t » ^\ - (4/2025)(Ax) 4 (the (Ax) 4 depen- 
dence was already observed in fl2]| ). In our simulations, 
Ax = 0.05, hence cui = 0.999998, approximately equal to 
2uj d . Therefore, the first peak in Fig. |] corresponds to 
the parametric resonance at half the value of the discrete- 
ness induced internal mode, very much like the resonance 
with the cf> 4 kink internal mode reported earlier B . 



B. Deformed initial kink 

In view of our negative result about the internal quasi- 
mode of the sG kink, we carried out more numerical sim- 
ulations specifically designed to find it. Following ||, we 
chose initial conditions given by 
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which correspond to a deformed kink if 1(0) =/= 1 or 
Z(0) 0. By starting their simulations with such ini- 
tial conditions, Boesch and Willis claimed that they were 
able to excite the internal quasi-mode, and hence we hope 
that, if it is indeed present, we should find it. 

At this point, we have to recall the work of Ma- 
jerm'kova, Gaididei and Braun who, as mentioned 
in the introduction, theoretically considered the prob- 
lem of a deformed kink given by the expressions above 
driven by a external constant force. When applied to this 
problem, our CC approach yields the same two ordinary 
differential equations as in Q for X{t) and l(t), [Eqs. 
@"(H) without the sine term, i.e., with only a constant 
force e]. We have solved these equations exactly for ar- 
bitrary initial conditions |l8| in terms of the Whittaker 
functions, concluding that the kink width oscillates only 
if Z(0) 7^ l s (the natural width of the undistorted kink, 
1/yl — u(0) 2 ) or 1(0) 7^ in (^|), whereas external force 
and dissipation (if present) only damp out those oscilla- 
tions. Hence, kink width oscillations should be easiest 
to detect without force and damping and, consequently, 
here we consider only free propagation of an initially 
distorted kink (a general analysis of the constant force 
problem will be presented elsewhere As a check of 

the validity of this procedure to probe internal modes, 
we have numerically verified that the </> 4 kink indeed be- 
haves in this fashion due to the excitation of the internal 
mode, exhibiting oscillations in the kink width very close 
to those predicted by the CC approach. We then turned 
to the sG case, and we found that the oscillations are not 
correctly described by the CC equations. In Fig. || we 
can see the evolution of the kink width computed from 
numerical simulation of the sG system, with ^(0) = 1, 
u(0) — 0, and 1(0) = 0.3. It is clear that l(t) is not a 
simple oscillatory function as in the <?i 4 case, and Fourier 
analysis yields the following frequencies: 014 = 0.9983, 
Qi = 1.0034, u) 2 = 1-0083, Q 3 = 1.0184, w 4 = 1.0335, 
w 5 = 1.0511, Cj 6 = 1.0712, oj 7 = 1.0963. As before, all 
these frequencies are very close to the discreteness in- 
duced internal mode and the first few radiational modes, 
with no evidence supporting the existence of the internal 
quasi-mode above the phonon band edge. 



IV. CONCLUSIONS 

The first conclusion of this work is that the existence 
of the proposed internal quasi-mode of the sG kink above 
the phonon band edge M is very unlikely. We have found 



neither resonances nor long-lived kink width oscillations 
arising from such a mode when the kink is subject to ac 
driving or to initial deformations, respectively. We are 
confident that we have carefully explored all the range 
of relevant frequencies, as we have accurately detected 
all possible resonances within that range. In addition, 
the comparison with the phenomenology observed for 
the cj) 4 kink, that possesses a true internal mode, shows 
that if the sG kink would have such an internal mode, 
we would have detected its influence through numeri- 
cal simulations as for the </> 4 equation In fact, 
now that we have seen that radiation modes are para- 
metrically excited by their corresponding half-frequency 
drivings, we believe we can suggest an explanation for 
the results of Boesch and Willis. Notice that they used 
Eq. (^) to calculate the lower phonon frequency, obtain- 
ing that wi = y/1 + (2tt/1000) 2 = 1.00002, whereas the 
internal "quasimode" frequency that they found in the 
numerical simulations was uj s = 1.004. As those two val- 
ues were largely different, they concluded that the quasi- 
mode was well above the lowest phonon frequency and 
different from it. However, as we have discussed in the 
previous sections, they used initial conditions in which 
the initial kink velocity was not zero; in this case, it can 
be easily shown that the corresponding frequencies of the 
phonons are given by 



ujk — k u(0) 
^1-^(0)' 



If we now insert the parameters used in [|| to perform a 
numerical experiment similar to the one discussed here, 
we find io\ — 1.004, i.e., the resonance observed by 
Boesch and Willis took place in fact with the lowest fre- 
quency phonon in the presence of a moving kink and not 
with any internal quasi-mode. 

The second conclusion of the present paper relates to 
the recently proposed internal modes induced by pertur- 
bations, in our case by discreteness [p^|-|l4|. By the two 
different types of numerical simulations used in our re- 
search, we have detected the internal mode bifurcating 
from the edge of the phonon band (of the continuum 
problem). The frequency of that mode agrees qualita- 
tively very well with the one predicted in [jl4j. Therefore, 
our research fully confirms the results of those papers 
and, in addition, shows that these perturbation induced 
internal modes behave very much like intrinsic ones. Our 
result, along with the finding of similar internal modes 
in other nonintegrable sG systems |19|, strongly supports 
the generality of the phenomenon predicted in |13| . Fi- 
nally, the fact that we are able to detect this novel mode 
reinforces our previous conclusion of the non-existence 
of the internal quasi-mode, as indeed the sG kink re- 
sponds to the ac driving or to deformations as it should 
if it had an internal mode, but for the discreteness in- 
duced one only. We thus believe that the present work 
definitely rules out the possibility of internal modes or 
"quasimodes" of sG kinks close to and above the lower 
phonon band edge. 



V. ACKNOWLEDGMENTS 



8.30 



We thank Yuri Kivshar and Yuri Gaididei for dis- 
cussions. Work at GISC (Leganes) has been sup- 
ported by DGESIC (Spain) grant PB96-0119. Travel be- 
tween Bayreuth and Madrid has been supported by "Ac- 
ciones Integradas Hispano-Alemanas" , a joint program of 
DAAD (Az. 314-AI) and DGESIC. 



kinter@math.uc3m.es 
anxo@math.uc3m.es 
franz.mertens@theo.phy.uni-bayreuth.de 
A. C. Scott, Nonlinear Science (Oxford University, Ox- 
ford, 1999). 

Yu. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 
763 (1989). 

A. Sanchez and A. R. Bishop, SIAM Rev. 40, 579 (1998). 
M. Peyrard and D. K. Campbell, Physica D 9, 33 (1983); 

D. K. Campbell, J. F. Schonfeld, and C. A. Wingate, ibid. 
9, 1 (1983); D. K. Campbell, M. Peyrard, and P. Sodano, 
ibid. 19, 165 (1986). 

Yu. S. Kivshar, F. Zhang, and L. Vazquez, Phys. Rev. 
Lett. 67, 1177 (1991); Phys. Rev. A 46, 5214 (1992). 
N. R. Quintero, A. Sanchez, and F. G. Mertens, Phys. 
Rev. Lett. 84, 871 (2000). 

S. Habib, G. Lythe, and F. G. Mertens, unpublished. 
R. Boesch and C. R. Willis, Phys. Rev. B 42, 2290 (1990). 
M. J. Rice, Phys. Rev. B 28, 3587 (1983). 
As a few examples, see P. Tchofo Dinda. R. Boesch, E. 
Coquet and C. R. Willis, Phys. Rev. B 46, 3311 (1992); 
Zhang Fei, Vladimir V. Konotop, Michel Peyrard and 
Luis Vazquez, Phys. Rev. E 48, 548 (1993); Fei Zhang, 
Phys. Rev. E 54, 4325 (1996); Physica D 110, 51 (1997). 

E. Majerm'kova, Yu. B. Gaididei, and O. M. Braun, Phys. 
Rev. E 52, 1241 (1995). 

O. M. Braun, Yu. S. Kivshar, and M. Peyrard, Phys. 
Rev. E 56, 6050 (1997). 

Yu. S. Kivshar, D. E. Pelinovsky, T. Cretegny, and M. 
Peyrard, Phys. Rev. Lett. 80, 5032 (1998). 
P. G. Kevrekidis and C. K. R. T. Jones, Phys. Rev. E 
61, 3114 (2000). 

F. G. Mertens, H. J. Schnitzer, and A. R. Bishop, Phys. 
Rev. B 56, 2510 (1997). 

N. R. Quintero and A. Sanchez, Eur. Phys. J. B 6, 133 
(1998). 

W. A. Strauss and L. Vazquez, J. Comput. Phys. 28, 271 
(1978). 

N. R. Quintero, A. Sanchez, and F. G. Mertens, unpub- 
lished. We note that in jll| only solutions for a special, 
restricted set of initial conditions were found. 
S. F. Mingaleev, Yu. B. Gaididei, E. Majerm'kova, and S. 
Shpyrko, Phys. Rev. E 61, 4454 (2000). 



8.25 



8.20 



8.15 



8.10 



8.05 



i.00 



ill 



0.45 0.48 0.50 



0.53 

5 



0.55 0.58 0.60 



FIG. 1. Mean energy of the sG system, in dimensionless 
units, when 8 is close to £Ir/2. Points have been calculated 
from numerical simulations for 10 000 < t < 25 000, with 
e = 0.01, 5 = tt/2, «(0) = 0, X(0) = 0. The parameters of 
the discretization are Ax = 0.05, At — 0.005 and L — 100. 
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FIG. 2. Evolution of the kink width, in dimensionless units, 
obtained from numerical simulations starting from a free, de- 
formed kink, with u(0) = 0, 1(0) = 1 and 1(0) = 0.3. The 
parameters of the discretization are L — 100, Ax = 0.1 and 
At = 0.01. 



Mode(n) 


w» 


LJ n /2 


0>n 


\cu n /2 - lj„\ 


1 


1.0019 


0.5009 


0.5009 


8.6 10~ 5 


2 


1.0079 


0.5039 


0.5040 


7 10~ b 


3 


1.0176 


0.5088 


0.5090 


1.9 10" 4 


4 


1.0311 


0.5155 


0.5161 


5.5 10~ 4 


5 


1.0482 


0.5241 


0.5250 


9.1 10~ 4 


6 


1.0687 


0.5343 


0.5356 


1.3 10~ d 


7 


1.0924 


0.5462 


0.5479 


1.7 10~ 3 


8 


1.1192 


0.5596 


0.5618 


2.2 10" d 



TABLE I. Comparison of the half of the first 8 frequencies 
<jj n /2 of radiational modes with the computed frequencies ui„ 
from numerical simulations of ([j]) at which peaks in E m arise 
(see Fig.0). 
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